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ABSTRACT 


The  nature  of  operations,  variety  of  missions  and 
conf igurations , and  requirement  for  rapid  response  of 
fleet  tactical  aircraft  has  given  rise  to  a unique  problem 
within  the  Naval  Aviation  community,  this  being  the  inter- 
face between  the  aircrews  and  the  aircraft  performance 
prediction  information  supplied  via  charts  and  data  in  the 
Naval  Air  Training  and  Operating  Procedures  Standardization 
(NATOPS)  Manuals.  The  data  presented  are  not  optimally 
used  due  to  the  method  of  presentation  and  the  time  and 
effort  required  to  extract  useful  predictions.  This  inves- 
tigation addressed  the  problem  by  developing  a suitable 
method  of  NATOPS  curve  presentation  that  has  sufficient 
simplicity  and  accuracy  for  application  to  current  compact 
computers.  The  procedure  is  defined  and  numerical  algor- 
ithsm  produced  which  demonstrate  the  feasibility  and  de- 
sirability of  this  type  of  implementation. 
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I . INTRODUCTION 

This  investigation  was  commenced  in  response  to  a 
generally  agreed  upon,  but  not  officially  defined,  fleet 
problem,  i.e.,  the  interface  between  operational  aircrews 
and  their  aircraft  performance  prediction  data.  Currently, 
in  the  rear  sections  of  all  fleet  aircraft  Naval  Air  Training 
Operations  and  Standardization  (NATOPS)  Manuals  [Ref.  5], 
an  extensive  set  of  performance  prediction  charts  are  included 
for  aircrew  use.  An  example  of  these  charts  for  the  A-7E 
aircraft  is  included  as  Appendix  A.  These  charts  represent 
performance  predictions  based  upon  interpolations  and  extrap- 
olations of  flight  test  data  compiled  at  the  Naval  Air  Test 
Center,  Patuxent  River,  Maryland.  They  are  included  in  the 
manuals  for  use  by  aircrews  in  the  prediction  of  aircraft 
performance  items  such  as  take-off  roll  and  velocities, 
cruise,  emergency  divert  (Bingo) , etc.  The  take-off  pre- 
diction parameters  starting  on  page  11-36  of  Appendix  A will 
be  used  as  a representative  example  and  have  been  prepared 
on  fold-out  sheets  for  ease  of  reference.  To  predict  take- 
off parameters,  the  aircrew  member  must  enter  the  chart  on 
page  11-36  of  the  NATCPS  Manual  with  runway  temperature  and 
pressure  altitude.  He  must  interpolate  as  necessary  between 
the  given  lines  and  proceed  as  in  the  dotted  line  example 
(refer  to  the  fold-out  of  Appendix  A numbered  11-36)  to 
the  chart  for  the  appropriate  double  datum  selection  (an 
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engine  parameter)  to  arrive  at  the  corresponding  take-off 
factor.  With  this  number  he  must  then  proceed  to  page  11-37 
and,  after  interpolating  for  gross  weight  correction,  pro- 
ceed as  indicated  to  arrive  at  the  basic  ground  roll  distance. 
Proceeding  to  page  11-38,  he  next  must  apply  a correction 
for  relative  humidity,  utilizing  the  relationship  in  the 
appropriate  paragraph,  followed  by  tracing  through  subsequent 
graphs  to  correct  for  relative  wind  (which  is  computed  from 
given  tower  winds  and  runway  heading),  runway  slope,  center 
of  gravity  position,  and  trailing  edge  flap  position.  Simi- 
lar graphical  procedures  exist  for  all  other  aircraft  per- 
formance predictions.  A compact  equivalence  of  these  charts, 
known  as  "Pocket  Check  Lists",  are  carried  by  all  aircrews 
in  flight  for  ready  reference.  In  actual  practice  individual 
squadrons  have  generated  numerous  NATOPS  data  cards  which 
represent  normally  encountered  conditions  and  were  derived 
from  the  NATOPS  charts.  In  the  take-off  data  prediction 
example  ten  input  parameters  were  required  to  compute  take- 
off roll  and  velocity.  Using  the  charts  has  proven  to  be 
cumbersome,  time  consuming,  and  prone  to  error.  The  use  of 
pocket  check  lists  is  more  convenient  but  less  accurate, 
and  the  individual  squadron  knee-board  cards  cover  only  a 
few  nominal  situations.  The  computation  of  take-off  roll 
is  relatively  simple  when  compared  with  determining  the 
parameters  to  establish  an  optimum  flight  profile  for  a 
complex  weapons  delivery  mission.  The  aforementioned  complex 
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problem  involves  a similar  use  of  cruise,  drag  count,  climb 

and  descent,  etc.  charts.  As  a further  illustration  of  the 

problem's  magnitude,  the  former  Director  of  the  School  of 

Aviation  Safety  at  the  Naval  Postgraduate  School,  Monterey, 

California,  has  repeatedly  required  that  each  member  of  a 

one-week  Command  Safety  Course  (prospective  squadron  commanding 

officers)  complete  the  weapons  delivery  problem  included 

here  as  Appendix  B.  Tnis  is  considered  as  a typical  fleet 

contingency  planning  problem,  requiring  that  in  one  hour 

these  experienced  aviators  compute  the  maximum  mission 

range  with  the  given  inputs  and  aircraft  parameters.  The 

results  of  these  tests,  given  to  a class  of  sixteen  such 

aviators,  were  described  by  the  Director  as  [Ref.  3] : 

"It  is  a startling  but  typical  fact  that 
the  correct  answer  of  538  NM  was  not 
achieved  by  any  member  of  the  class. 

The  closest  answer  was  in  error  by  126 
miles  and  the  spread  of  answers  ranged 
from  336  to  868  NM.  Additionally  the 
correct  answer  was  attained  by  the  class 
instructor  only  after  a measured  16  hours 
of  effort  with  the  NATOPS  Manual." 

The  preceding  example  serves  as  a graphic  illustration  of 

a widely  accepted  problem,  the  elements  of  which  may  be 

summarized  as: 

1.  The  only  complete  and  accurate  performance 
prediction  data  available  to  fleet  aircrews  are 
those  provided  by  the  performance  chart  sections 
of  the  appropriate  aircraft  NATOPS  Manuals. 
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2.  Experienced  fleet  aviators  do  not  use  the 
source  information  for  daily  computation  of 
NATOPS  problems  with  enough  regularity  to 
insure  competency. 

3.  Personal  experience  and  interviews  with 
fleet  aviators  indicated  that  the  NATOPS 
performance  presentation  contained  highly 
useful  and  vital  data  but  was  too  cumbersome, 
and  data  extraction  too  time  consuming,  for 
regular  operational  use. 

The  purpose  of  the  studies  described  herein  was  to 
examine  an  alternate  approach  for  assisting  the  aircrew 
member  to  solve  the  various  aircraft  performance  problems 
which  are  required  for  safe  and  effective  operation  of  a 
Naval  aircraft  in  fulfilling  its'  mission.  The  subsequent 
discussion  will  proceed  with  the  assumptions  made,  the 
method  used  for  the  solution,  and  the  results  of  this  study. 
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II. 


DISCUSSION 


A.  GUIDELINES 

At  the  onset  of  the  investigation  a set  of  guidelines 
was  adopted  to  limit  its  scope  and  provide  a definitive 
goal.  This  was  deemed  necessary  because  of  the  controversy 
within  the  aviation  community  over  all  facets  of  NATOPS 
including  the  validity  of  some  of  the  data.  It  was  beyond 
the  scope  of  this  investigation  to  question  that  validity 
or  investigate  any  reported  discrepancies . ^ The  most 
current  NATOPS  data  available  are  the  official  fleet  stan- 
dard and  it  was  merely  the  representation  of  those  data 
that  inspired  this  study.  Thus,  from  the  above  discussion 
the  first  guideline  was  defined: 

1.  Any  method  of  representation  of  NATOPS  data  must 
be  compared  to  the  officially  published  NATOPS  data 
for  accuracy  and  completeness. 

Throughout  the  life  cycle  of  an  aircraft,  numerous 
changes  to  the  NATOPS  procedures  and  performance  data  are 
published  which  represent  revisions  due  to  lessons  learned 
through  extensive  operating  and  testing.  The  realization 
of  this  fact  led  to  guideline  number  two: 


During  the  course  of  this  investigation  several  such 
discrepancies  were  found  and  changes  submitted  in  accordance 
with  existing  change  submittal  procedures. 
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2.  To  remain  useful  and  current,  any  system  of 
data  presentation  must  be  adaptable  to  changes  as 
they  are  promulgated. 

The  mission  variety  and  uniqueness  of  carrier  based, 
single-seat  aircraft,  especially  under  emission  control 
conditions,  make  it  highly  desirable  that  the  information 
contained  in  NATOPS  be  readily  available  to  the  pilots  while 
in  the  aircraft.  Airborne  workload,  lack  of  additional 
cockpit  space,  and  lack  of  time  in  critical  situations  led 
to  the  third  guideline: 

3.  The  data  presentation  system  utilized  must  be 
light,  small,  simple  to  use,  reliable,  and  require 
a minimum  amount  of  pilot's  attention  during  use 
to  insure  acceptability  from  the  standpoints  of 
mission  accomplishment  and  flight  safety. 

Finally  a standard  of  measurement  was  defined  for 
accuracy  comparison  with  NATOPS  data.  These  figures  were 
chosen  based  upon  tha  experience  of  the  author,  and  were 
used  to  determine  when  the  degree  of  complexity  of  the 
computer  solution  was  adequate  for  this  application. 
Guideline  four  lists  the  desired  computational  tolerances. 

4.  One  knot  of  airspeed 

100  feet  of  altitude  of  ground-roll  distance 
100  pounds  of  weight 
Ten  seconds  of  time 
100  pounds  of  fuel 
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After  defining  the  above,  an  investigation  was  commenced 
toward  utilizing  digital  computers  as  a solution,  with  the 
goal  being  the  reduction  of  the  data  contained  in  NATOPS 
into  an  algorithm  for  use  in  desk  , "hand-held",  and/or 
existing  on-board  aircraft  calculators  and  computers. 
Informal  liaison  with  senior  and  junior  attack  and  fighter 
fleet  aviators  has  indicated  strong  support  for  a solution 
of  this  type  when  restrained  to  the  stated  guidelines. 

B.  COMPUTATIONAL  METHOD  USED 

Procedurally  the  entire  project  was  subdivided  as 
follows : 

1.  Investigation  of  a numerical  technique  for  general 
computer  use. 

2.  Application  of  thst  numerical  technique  to 
specific  aircraft  performance  data 

3.  Investigation  of  hardware  and  implementation  of 
the  results  of  part  2,  above,  for  fleet  use. 

This  investigation  was  principally  concerned  with  step  1, 
above,  while  subsequent  and  concurrent  investigations  were 
proceeding  with  the  work  in  the  other  areas.  Two  basic 
techniques  were  investigated: 

1.  Data  Storage  t 

2.  Computational  technique 
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1 

I 

Data  storage,  while  mathematically  trivial,  was 
, determined  to  be  unacceptable  for  this  application  due  to 

the  amount  required  to  insure  completeness.  An  examination 
of  a set  of  curves  (e.g.,  take-off)  showed  a near  infinite 

|i 

number  of  solutions  for  the  ten  input  variables  used.  The 
current  squadron  method  for  bypassing  NATOPS  by  using  knee- 
board  cards  for  quick  reference  to  nominal  solutions  corres- 
ponds to  the  data  storage  method  on  a much  smaller  scale. 

While  many  more  data  poinrs  could  be  stored  utilizing  calcu- 
lators, it  was  felt  that  a table  lock-up  solution  would  lead 
to  incorrect  interpolation  and  extrapolation  when  seeking 
information  utilizing  parameters  other  than  the  ones  stored. 
Interpolation  routines  could  be  utilized  when  this  occurred, 
but  each  routine  would  vary  depending  upon  the  relationship 
of  the  data  points.  In  many  cases  the  data  vary  in  a differ- 
ent non-linear  manner,  thus  causing  the  number  of  routines 
used  to  be  unacceptably  large  in  order  to  produce  the  desired 
results.  Additionally,  the  number  of  storable  points  is 
limited  by  computer  size.  Fleet  tactical  aircraft  continually 
operate  with  widely  varying  fuel,  gross  weight  and  drag 
coefficient  values.  These  parameters  greatly  affect  nearly 
all  computations  made.  The  data  storage  method  for  mission 
planning  purposes  would  require  a near  infinite  number  of 
precomputed  solutions  using  these  parameters,  and  would 
therefore  be  too  complex  for  use.  Finally  the  computer 
state-of-the-art  was  sufficiently  advanced  to  and  even 
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beyond  the  capability  for  computing  the  performance  of 
any  number  of  inputs,  utilizing  a fraction  of  the  memory 
space  required  for  data  storage. 

The  data  storage  method,  after  a brief  trial,  was 
not  investigated  further.  A summary  of  the  reasons  for 
this  decision  would  include: 

1.  The  number  of  variables  involved  would  require 
too  much  data  storage  space  to  produce  the  desired 
result. 

2.  A reduced  set  of  data  would  be  of  little 
value  in  mission  planning  situations  due  to  the 
accuracy  of  the  required  results. 

3.  Inaccurate  (and  thus  unsafe)  interpolation  and 
extrapolation  of  the  stored  data  would  be  a 
natural  out-growth  of  this  type  of  system 
(similar  to  the  current  use  of  knee-board  cards) . 

4.  Computer  programming  state-of-the-art  was 
far  advanced  beyond  the  requirement  for  this 
type  of  method. 

The  computational  method  was  considered  to  be  well 
within  the  capability  of  "hand-held"  calculators  from  both 
memory  size  and  computational  standpoints.  The  basic 
steps  for  the  implementation  of  this  method  are  listed 
below : 
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1.  Reduce  the  individual  curves  or  sets  of 
data  to  analytical  expressions. 

2.  Coordinate  these  expressions,  of  each  of 
the  above,  to  output  the  dependent  variable 
for  each  chart  or  family  of  curves,  providing 
for  interpolation  between  the  curves. 

3.  Coordinate  the  inputs  and  individual  curve 
outputs  into  a computer  program  for  validation 
and  pilot  use. 

The  method  used  for  reducing  the  curves  and  charts  to 
analytical  expressions  was  simply  the  Least  Squares  Fit 
method.  A precise  explanation  of  this  method  can  be  found 
in  Refs.  1 and  2.  Appendix  C contains  a brief  explanation 
and  sample  application  of  the  technique.  The  majority  of 
the  curves  or  data  sets  found  in  NATOPS  could  be  sufficiently 
described  by  first,  second,  or  third  degree  polynomials. 

This  fact  was  assumed  from  intuition  and  supported  by  numer- 
our  iterative  curve  fits.  A variety  of  computer  programs 
were  in  existence  to  aid  in  producing  the  desired  solutions. 
The  Hewlett-Packard  (HP)  9330  desk  calculator  had  available 
one  such  program  and  was  used  for  this  analysis  (c.f., 
Reference  6 for  a detailed  description  and  set  of  instruc- 
tions for  use  of  this  technique) . Prior  to  commencement  of 
the  data  reduction,  photographically  enlarged  versions  of 
the  charts  under  investigation  were  obtained  to  aid  in  the 
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precise  recording  of  information.  Using  these  charts  the 
following  procedure  was  utilized  to  produce  the  analytical 
expressions  for  each  individual  curve: 

1.  Record  data  points  from  each  curve  of 
a family  of  curves  presented  on  the  chart 
(Refer  to  the  foldout  charts  for  take-off 
prediction  provided  in  Appendix  A) . 

2.  Utilize  these  points  in  the  HP  polynomial 

\ 

regression  program  to  generate  the  coefficients 
of  the  analytical  expression  of  sufficient 
degree  for  the  curve  being  analyzed. 

3.  Use  the  coefficients  generated  to  check 
polynomial  accuracy  when  compared  to  the 
NATOPS  curves. 

4.  Iterate  the  procedure,  reducing  the  degree 
of  the  polynomial  in  an  attempt  to  simplify  the 
expressions  and  retain  sufficient  accuracy  to 
comply  with  the  specifications  of  guideline  4. 

Once  a set  of  equations  was  obtained  for  a family  of 
curves,  an  investigation  was  commenced  to  determine  a 
simple  method  of  coordinating  the  curves,  and  providing 
a method  of  producing  an  output  for  those  input  parameters 
which  might  fall  between  the  curves  of  a family.  This  pre- 
sented the  greatest  problem.  The  curves  of  the  top  chart 
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on  page  11-38  of  Appendix  A were  used  as  an  example  (fold- 
out  provided) . Given  the  set  of  five  polynomial  expressions 
corresponding  to  the  respective  curves  for  runway  pressure 
altitude,  a linear  interpolation  routine  was  produced  and 
tested  to  predict  values  between  the  curves.  This  routine 
represented  the  first  attempt  at  providing  a means  of 
accurately  computing  values  between  the  curves.  The  inde- 
pendent variables  were  the  pressure  altitude  and  the  runway 
temperature.  Using  the  interpolation  routine  for  a pressure 
altitude  of  2,595  feet  and  95  °F,  for  example,  the  program 
first  calculated  the  output  from  the  2,000  foot  curve  poly- 
nomial. It  then  did  the  same  for  the  4,000  foot  polynomial 
and  linearly  interpolated  to  find  a value  for  the  original 
input  of  2,595  feet.  This  procedure  produced  unacceptable 
results  for  the  reasons  stated  in  the  following: 

1.  The  linear  interpolation  was  not  a true 
model  of  the  actual  variance  between  the  curves. 

A visual  inspection  revealed  that  the  spacing 
between  the  curves  increased  with  altitude.  An 
attempt  was  made  to  correct  the  effect  by  multi- 
plying the  result  by  a complicated  bias  factor  with 
analytically  acceptable  results.  The  method  used 
to  arrive  at  this  factor  was  highly  iterative  and 
time  consuming  and  not  considered  an  acceptable 
procedure  considering  the  number  of  curves  to 
be  analyzed. 
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2.  Additionally  the  program  required  for  selecting 
the  curves  for  interpolation  for  a given  input, 
complete  with  interpolation  routine,  was  considered 
too  lengthy  and  repetitive  for  small  calculator 
use.  This  would  be  especially  critical  if  the 
family  of  curves  contained  many  individual  curves, 
as  do  those  on  page  11-38  of  Appendix  A. 


Given  that  the  method  of  linear  interpolation  between 
curves  of  a family  was  unacceptable  for  the  reasons  above, 
an  alternate  method  of  expressing  a family  of  curves  analy- 
tically was  investigated.  This  method  was  suggested  by 
Lt.  Keith  Winters  [Ref.  4].  It  will  be  explained  using 
the  example  contained  in  the  following  paragraphs. 

The  problem  was  to  take  a set  of  polynomials  of  the 
same  degree  representing  a family  of  curves  and  reduce  these 
to  an  equation  or  small  set  of  equations  that  would  generate 
a solution  for  these  polynomials  or  any  value  in  between 
the  curves  they  represent. 

Given  was  a set  of  equations  for  the  five  curves  of  the 
previous  example  (all  of  the  same  degree)  of  the  form 


Y1  A11  + A12X  + A13X  A1NX 


Y2  = A21  + A22X  + A23x2  A2NxN_1 


’ 2 N-l 

Y = A , + A 0X  + A ,X  A X 

m ml  m2  m3  mN 
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or  more  specifically  for  this  example 

Y1  = A11  + A12X  + A13x2  + A14x3  0 feet 

pressure 

altitude 

Y 2 = A21  + A22X  + A23X2  + A24X3  1000  feet 

Y3  = A31  + A32X  + A33/2  + A34x3  2000  feet 

Y4  = A41  + A42X  + A43x2  + A44X3  4000  feet 

Y-  = Ac  + Ac-X  + AcoX2  + A- • X3  8000  feet 

5 5;  52  53  54 

The  A^  terms  were  the  polynomial  coefficients  produced 
by  the  Least  Squares  Fit  program  for  each  curve  of  chart  1. 
The  Y terms  were  the  fictitious  numbers  assigned  as  the 
output  dependent  variables  of  chart  1 (input  to  chart  2) 
and  X was  the  independent  variable,  temperature,  of  the 
same  example.  Again  it  was  desired  that,  given  any  altitude 
and  temperature  within  the  chart  limits,  a Y could  be  com- 
puted precisely.  To  do  this  the  program  must  first  take  the 
given  altitude  and  generate  either  computationally,  or  from 
memory,  the  coefficients  of  the  corresponding  third  degree 
polynomial.  To  provide  for  the  computation  of  these  coeffi- 
cients separate  plots  were  made  of  the  A ^ , A^2 , Am3  and 
A^^  coefficients  using  their  corresponding  input  altitude 
values  as  abscissae.  Figure  1 illustrates  the  concept. 
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Altitude  (Feet) 


Figure  1.  Altitude  versus  polynomial  coefficients 
for  a series  of  third  degree  equations 

Since  each  curve  contained  five  data  points,  a fourth 
degree  (or  co-location)  polynomial  was  plotted  and  the 
coefficients  produced  by  the  HP  9830  program.  With 
these  coefficients  (Bed)  equations  of  the  form  below 
were  produced: 
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These  equations  when  used  with  the  following  equation 


Y = A ■,  + A -X  + A ,X2  + a X3 
ml  m2  m3  m4 


would  produce  the  desired  output  for  any  value  of  temperature 
and  altitude.  The  following  chart  of  symbol  definitions 
was  included  to  ensure  clarity: 


M = Any  value  of  pressure  altitude  between  0 
and  8000  feet 

x = Any  value  of  temperature  between  0 and  120 °F 

Y = The  fictitious  number  representing  the  output 

of  chart  1 for  any  given  input 

Bed  = Coefficients  of  the  polynomials  used  to  predict 

the  corresponding  coefficient  of  the  Y polynomial 
for  any  given  input  pressure  altitude  (M) 

Amn  = The  coefficient  of  the  x polynomial  used  to 
generate  the  output  number  of  chart  1 from 
the  given  temperature. 


The  task  of  the  person  attempting  to  program  a set  of 
NATOPS  curves  using  this  technique  is  summarized  below: 
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1.  Generate  the  coefficients  for  each  curve  by: 

a.  Recording  data  points  from  the  curves. 

b.  Utilize  a computer  solution  or  manually 
fit  the  points  to  output  a polynomial  of 
the  same  degree  representing  each  curve. 

2.  Plot  the  similar  coefficients  of  each  curve 
as  a dependent  variable  against  the  corres- 
ponding independent  variable  as  was  done  in 
Figure  1. 

3.  Obtain  polynomial  expressions  of  sufficient 
degree  for  each  of  these  coefficient  curves. 

4.  Write  a computer  algorithm  that  will  perform 
the  following: 

a.  Utilize  the  polynomials  of  step  3 to 
generate  the  coefficients  of  the  polynomials 
of  step  2. 

b.  Utilize  the  coefficients  generated  in  step 
4a.  in  a polynomial  that  will  output  the 
independent  variable  of  a family  of  curves. 

These  steps  are  a general  description  of  the  successful 
procedure  which  was  generated  by  this  investigation.  A 
detailed  example  using  actual  numbers  has  been  included  in 
the  next  section  for  reference  purposes.  Throughout  this 
explanation,  emphasis  has  been  placed  upon  the  need  for 
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4 


interpolating  accurately  for  any  value  within  the  given 
curves.  The  method  used  also  provides  for  mathematic 
extrapolation  outside  the  limits  of  the  given  NATOPS  curves. 

Care  was  taKen  not  to  provide  for  this  extrapolation  in 
any  of  the  programs  generated  by  this  technique.  This  was 
done  either  by  limiting  the  range  of  the  independent  varia- 
bles which  the  program  would  accept  to  values  within  the 
stated  NATOPS  limits,  or  by  ensuring  that  the  output  of 
each  separate  calculation  did  not  exceed  the  corresponding 
chart  limits.  Mathematical  extrapolation  of  NATOPS  data, 
without  flight  test  verification  and  approval  from  the  Naval 

I 

Air  Test  Center,  can  lead  to  a computer  prediction  that  does 

not  correspond  to  aircraft  performance.  While  precise  \ 

interpolation  is  well  within  the  capability  of  the  process 
developed  by  this  investigation,  extrapolation  in  this  case 

> 

is  dangerous  and  was  not  within  the  intent  of  the  study. 

C.  DETAILED  EXAMPLE  OF  THE  COMPUTATIONAL  METHOD 

The  computational  steps  listed  on  page  24  were  performed 
for  the  take-off  charts  of  the  A7E  NATOPS  manual  and  were 
included  as  an  example.  The  number  designation  of  each 
step  corresponds  to  the  instruction  as  listed.  The  problem 
was  to  generate  a set  of  equations  that  would  describe  the 
family  of  curves  on  page  11-36  (Top  chart)  of  Appendix  A 
and  provide  for  interpolation  between  all  curves  to  output 
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the  dependent  variable  to  be  used  in  the  subsequent  chart 
on  that  page. 

la.  The  data  shown  below  were  extracted  from 
the  top  chart  of  page  11-36  Appendix  A. 
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lb.  These  data  were  entered  into  the  afore- 
mentioned HP  9830  computer  program, 
producing  the  following  coefficients  of 
a third  degree  polynomial  of  the  form 


Y = A 


ml 


+ A _X 
m2 


+ A ,X2  + A .X  3 
m3  m4 


where 

m = altitude 
X = temperature 


2 

It  may  be  noted  that  the  chart  used  does  not  have  printed 
values  for  the  dependent  variable.  In  NATOPS  curves  these  are  not 
required  since  the  normal  procedure  is  to  proceed  graphically  from 
chart  to  chart.  For  the  numerical  analysis  of  this  study  a set  of 
fictitious  numbers  was  assigned  to  the  dependent  variable  of  this 
chart  and  to  the  independent  variable  of  chart  2 page  11-36. 
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Coefficients 

> 

3 
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o 
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- MX« 

2.  The  coefficients  were  then  plotted 
against  altitude  as  depicted  in 
Figures  2,  3,  4,  and  5. 


8c?oo 


> 

Altitude 
Figure  4:  versus  Altitude 


Altitude 


Figure  5: 


Am4  versus  Altitude 


3.  Using  the  HP  9830  computer  program  as 

above  the  coefficients  of  each  of  these 
third  degree  curves  were  generated  and 
placed  in  the  following  polynomials 
of  the  form: 

2 3 

= B . + B _m  + B ..m  + B .m 


where 

m = Altitude  (Feet) 

A . = 13.1  - 1.7x10  4m  - 2.1x10  7m^  + 3.7x10 
mi 

Am2  = 4 . 5x1 0-2  - 7.9xl0'6 a + 3.8xl0'9m2  - 9.7xl0_13m3 

Am3  = 1 . 3xlO~3  - 3.2xl0_7m  + 4.1xl0~10m2  - 3.5xl0_14m3 

A^  = 1.9x10  3 + 1.4x10  ^m  - 9.5x10  ^2m2  + 2.0x10  33m3 

4a. b.  The  enclosed  computer  program  (following 
the  Appendix)  satisfies  the  requirements  of 
steps  4. a and  b,  page  24.  It  coordinates 
all  the  input  variables  needed  for  computing 
take-off  velocity  and  ground  roll  distance, 
and  outputs  those  values  within  the  accuracy 
required  by  Guideline  4 of  page  13. 
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III. 


RESULTS 


The  results  of  this  investigation  were  that: 

1.  A technique  for  reducing  families  of 
curves  found  in  NATOPS  manuals  to  simple 
equations  that  provide  for  proper  interpolation 
was  developed  and  tested. 

2.  A sample  program  for  an  entire  performance 
example  was  prepared  to  illustrate  the  use  of 
the  polynomials  in  generating  a solution  to  a 
typical  performance  problem. 

The  computational  technique  was  utilized  in  providing 
coefficients  for  polynomials  for  all  families  of  curves 
required  to  compute  take-off  roll  and  velocity  for  the  A7E 
aircraft.  A list  of  the  program  steps  was  included  as  an 
enclosure  to  this  thesis.  The  curves  of  the  charts  for 
the  ten  input  parameters  were  treated  in  a manner  similar 
to  the  example  of  the  previous  section. 

As  can  be  seen  from  the  equations  of  each  section  of  the 
enclosed  program  the  number  of  polynomial  expressions  re- 
quired to  describe  a family  of  curves  is  not  related  to  the 
number  of  curves  within  the  family.  Rather,  it  is  a function 
of  the  degree  of  the  curves  within  the  family.  If  the  degree 
of  the  curves  were  represented  by  the  letter  w,  then  the 
number  of  equations  required  (s)  would  be 
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s = w + 2 

This  fact  caused  significant  reduction  in  the  size,  and 
an  increase  in  the  computational  accuracy,  of  the  final 
program.  To  illustrate  an  extreme  example,  the  family  of 
curves  associated  with  computing  the  take-off  distance  wind 
correction  will  be  referenced.  It  can  be  seen  from  this 
family  {Located  in  Appendix  A,  page  11-38)  that  fifteen 
second  degree  curves  were  utilized  to  describe  this  correc- 
tion. The  total  number  of  curves  required  to  compute  the 
correction  using  the  method  of  this  investigation  was 

s = 2 + 2 = 4 . 

This  included  the  three  equations  for  computing  the  coeffi- 
cients and  one  equation  for  coordinating  the  coefficients 
and  computing  the  dependent  variable  from  the  output  of 
the  previous  algorithm.  It  is  worthwhile  at  this  point  to 
emphasize  that  these  four  equations  were  not  only  valid  for 
the  data  that  formed  the  given  curves,  but  also  for  any 
intermediate  values  retaining  the  same  accuracy.  Alternate 
methods  of  using  equations  for  all  fifteen  curves,  plus  a 
discrimination  and  interpolation  routine,  would  produce 
less  accurate  results  and  require  more  computer  space. 

The  sample  program  for  take-off  velocity  and  distance 
computation  was  implemented  for  demonstration  purposes  for 
the  HP  9830.  The  program  and  computer  combination  met  all 
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the  requirements  set  forth  in  the  guidelines  with  respect 
to  ease  of  use,  accuracy,  speed,  reduced  memory  required, 
and  all  around  desirability.  Informal  tests  with  fleet 
aviators  demonstrated  that  the  program  produced  the  same 
results  as  the  charts  but  required  less  time  by  an  order 
of  magnitude,  and  little  or  no  user  effort.  Implementation 
of  the  entire  set  of  NATOPS  curves  on  the  HP  9830  desk 
computer  would  be  unsuitable  for  aircraft  use  but  could 

provide  preflight  and  contingency  planners  with  an  invaluable  , 

1 

tool  for  ground  and  shipboard  use. 

Concurrent  investigations  were  performed  in  this  area 
by  Major  J.D.  Restivo,  USMC,  and  Lt  G.L.  Roger,  USN,  at  the 
Naval  Postgraduate  School.  Major  Restivo  applied  the  tech- 
nique of  this  thesis  to  the  A6  aircraft  NATOPS  curves,  while 
Lt  Roger  implemented  the  entire  set  of  A7E  curves  in  a 
"hand-held"  calculator,  namely,  the  Texas  Instruments 
TI  59. 

With  the  completion  of  this  thesis,  the  feasibility 
of  reducing  the  NATOPS  performance  data  and  curves  has 
been  demonstrated,  and  a procedure  documented  for  the 
process.  The  responsibility  for  further  investigation  in 
this  area  has  been  assumed  by  Commander,  Light  Attack  Wings, 

Pacific  (COMLATWINGSPAC) , Naval  Air  Station,  Lemoore, 

California.  To  this  end  the  following  conclusions  and 
recommendations  are  included  as  the  final  section. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  is  concluded  from  this  investigation  that  adaptation 
of  the  NATOPS  performance  data  to  digital  computers  for 
aircrew  use  was  entirely  feasible.  Furthermore,  a simple 
procedural  method  was  found  and  demonstrated  for  reducing 
a large  family  of  curves  to  a minimum  number  of  polynomial 
expressions  which  could  be  coordinated  in  a program  for 
precisely  computing  aircraft  performance  for  any  variation 
of  the  input  parameters.  Finally,  it  was  shown  that  the 
computerized  method  of  NATOPS  presentation  produced  a quick 
and  easy  method  for  aircrews  and  planners  to  predict 
aircraft  performance  for  missions  and  contingencies. 

The  following  recommendations  are  made  as  a result  of 
this  study: 

1.  The  activity  which  has  been  assigned  the 
task  of  managing  the  research,  development, 
and  procurement  for  use  by  tactical  aircraft 
crews  should  consider  the  procedures  developed 
in  this  thesis  study  for  achieving  the  end 
result . 

2.  Funds  should  be  provided  for  a coordinated 
Navy-wide  effort  to  ensure  standardization  of 
technique  and  hardware  to  reduce  redundant 
efforts  in  this  area  of  rapidly  increasing 
interest. 
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NAVAIR  01-45AAE-1 


APPENDIX  A 


IAKEOFF  FACTOR  (A-7E) 


ENGINE  TF41  A 2 
FUEL  GRADE  JP-5 
FUEL  DENSITY  6.8  LB/GAL. 


MODEL  A-7E 

DATA  BASIS  FLIGHT  TEST 
DATE  NOVEMBER  1971 


JbFUt.  -04-  74 


NAVAIR  01  45AAE1 


TAKEOFF  GROUND  ROLL  DISTANCE  (A-7E) 


MODEL:  A-7E 

DATA  BASIS  FLIGHT  TEST 
DATE:  NOVEMBER  1971 


CONDITIONS: 

LEVEL  HARD  SURFACE  RUNWAY 
MILITARY  RATED  THRUST 
LANDING  CONFIGURATION 
ZERO  HEADWIND 
CG:  26%  MAC 
FULL  FLAPS 


ENGINE:  TF41-A-2 
FUEL  GRADE:  JP-5 
FUEL  DENSITY:  6.8  LB/GAL. 


1 


tfNOTC* 


For  minimum  ground  roll  corresponding  to  minimum 
lift-off  speed,  subtract  500  feet. 


For  humidity  effects  on  takeoff  distance,  ground  roll 
distances  should  be  increased  1%  for  each  10%  increase 
in  the  relative  humidity  above  40%. 


7»»F?H  /(  t ) -04  - / 4 


Figure  11-12  ( Sheer  1) 
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NAVAIR  01  45AAE-1 
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TAKEOFF  GROUND  ROLL  DISTANCE  (A-7E) 


ADJUSTED  GROUND  ROLL 

MODEL  A 7E 

DATA  BASIS.  FLIGHT  TEST 
DATE  NOVEMBER  1971 


DISTANCE 

CONDITIONS: 

HARO  SURFACE  RUNWAY 
MILITARY  RATED  THRUST 
LANDING  CONFIGURATION 
LEADING  EDGE  FLAPS  DOWN 


ENGINE  TF41-A-2 

FUEL  GRADE  JP-5 

FUEL  DENSITY:  6.8  LB /GAL 


* NOTE* 


For  humidity  effects  on  takeoff  distance,  ground  roll 
distances  should  be  increased  1%  for  each  10%  increase 
in  the  relative  humidity  above  40%. 
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Figure  11-12  (Sheet  2) 
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APPENDIX  B 


A- 7 NATOPS  PROBLEM 


An  international  incident  has  occurred  which  necessi- 
tates the  rapid  protection  and  extraction  of  U.S.  military 
personnel  within  the  strife  stricken  country. 

You  are  the  Commande  r>f  an  A-7E  squadron  and  have  been 
given  the  task  of  planning  and  conducting  the  missions  in 
support  of  the  extraction  of  the  U.S.  personnel.  Given  the 
following  information,  determine  the  greatest  distance 
from  the  target  at  which  you  could  launch  and  return  to  the 
same  take-off  point. 

1.  Fuel  Load  - 14,000#  (2,000#  in  each  aero  I D) 

2.  Ordnance  Load 

a.  1 AGM-45  (Skrike)  on  Sta  #1  and  Sta  #8 

b.  3 MK-20  (Rockeye)  on  Ter  on  Sta  #2 

c.  1 Aero  I D (Full  JP5)  on  Sta  #3  and  Sta  #6 

d.  1 AIM  9D  (Sidewinder)  on  Sta  #4  and  Sta  #5 

e.  3 MK-83  ( LDGP ) on  Ter  on  Sta  #7 

3.  Basic  Aircraft  Weight  - 21,000# 

4.  Area  temperature  at  sea  level  - 30°C 

5.  Use  300#  fuel  for  Start/Taxi/Catapult 

6.  Rendezvous  at  10,000'  and  25  NM  - use  800#  fuel  during 
rendezvous 

7.  Descend  from  cruise  altitude  to  target 

a.  Descend  to  12,000' 

b.  Use  2,700#/Hr  fuel  flow 

c.  Maintain  370  Kts  ground  speed 

d.  Begin  descent  at  40  NM 

8.  Enroute  winds  are  forecast 

a.  To  the  target  - 40  Kts  headwind 

b.  Return  to  ship  - 70  Kts  tailwind 

c.  Descents  - no  wind 
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9. 


Use  2,000#  fuel  in  target  area 


10.  Expend  the  following  ordnance 

a . 2 AGM- 4 5 

b.  3 MK-20 

c.  3 MK-83 

d.  Retain  all  other  racks,  ordnance,  tanks 

11.  Climb  to  cruise  altitude  for  return  will  be  from  sea 
level. 

12.  Plan  for  Charlie  on  arrival  with  1,000#  internal  fuel 
on  the  ball. 

The  following  charts  and  tables  have  been  taken  from  the 
NATOPS  for  your  use. 
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APPENDIX  C 


Least  Squares  Fit  Approximation 

References  1 and  2 describe  the  Least  Squares  Fit 
Approximation  in  detail.  In  general  the  problem  is  to 
represent  a set  of  "n"  data  points  in  two  dimensional  space 

X . , Y . i = 1 to  n 
l l 

\ 

by  a polynomial  expression  of  a curve  whose  degree  is  less 

than  "n".  Two  classes  of  problems  exist:  ; 

1.  Linearly  independent  - Those  in  which  the 
degree  (d)  of  the  polynomial  is  one  less 
than  the  number  of  data  points 

d = n - 1 (1)  j 

2.  Linearly  dependent  - Those  in  which  the  degree 
(d)  is  less  than  n-1 

d < n - 1 (2) 

As  an  example,  a set  of  four  (4)  data  points  randomly 
spaced  was  chosen.  If  a third  degree  polynomial  of  the 
form 


l 
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Y = A + BX  + CX2  + DX3  (3) 

were  desired,  and  the  data  points  X^  and  Y ^ were  inserted 
(i  = 1 to  4)  into  four  such  equations,  an  exact  solution 
for  the  four  unknown  coefficients  would  exist.  These  four 
unknowns  could  be  found  from  the  four  equations  by  numerous 
conventional  techniques  (Direct  substitution,  Cramer's  rule, 
etc. ) . The  polynomial  expression  generated  would  be  termed 
a " col- location"  polynomial  because  its  plot  would  pass 
through  all  data  points. 

It  is  often  advantageous  to  describe  a set  of  data  points 
by  a curve  that  does  not  pass  through  each  point.  This 
type  of  polynomial  would  be  termed  a "regression"  equation. 

For  any  set  of  data  points  an  infinite  number  of  regression 
expressions  exist  for  any  specified  degree  (except  the 
linearly  independent  case)  and  the  object  of  the  Least  Squares 
Method  is  to  find  the  polynomial  coefficients  of  the  chosen 
degree  that  best  describe  the  data  points.  In  the  previous 
example  of  four  data  points,  assume  that,  instead  of  the 
third  degree  form  chosen,  a second  degree  equation  were 
selected  of  the  form 

Y = A + BX  + CX2  (4) 

With  four  data  points,  the  polynomial  is  overspecified  and 
thus  linearly  dependent.  For  this  case  an  infinite  number 
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of  solutions  exist  for  the  coefficients  a,  b and  c.  If 
an  error  term  (6)  were  defined  for  any  given  X,Y  pair  as 

$1  = |YL  - A + + CXl2|  (5) 

a total  squared  error  term  (E)  could  then  be  defined  by 
squaring  and  summing  the  terms  attained: 

N 

E = l 6X2  (6) 

i=l 

If  E were  then  minimized  for  any  given  degree  chosen,  the 
best  Least  Squares  Fit  would  have  been  achieved. 

If  the  values  for  6 from  equation  5 were  inserted  in 
equation  6 and  the  partial  derivative  of  E were  taken  with 
respect  to  the  coefficient  A an  equation  would  be  generated 
that  when  set  equal  to  zero  (0)  would  define  a minimum 
value  of  E for  a given  value  of  A.  If  the  same  operation 
were  performed  with  respect  to' the  coefficients  B and  C then 
three  equations  would  be  generated  with  three  unknowns  (A, 

B,  and  C) . The  solution  of  these  simultaneous  equations 
would  produce  the  coefficients  A,  B,  and  C,  that  would 
minimize  the  value  of  E and  hence  would  produce  a Least 
Squares  Fit  approximation  to  a set  of  linearly  dependent 
equations . 

A numberical  procedure  has  been  developed  to  accomplish 
this  task.  An  example  of  this  procedure  has  been  included 
in  the  following  paragraphs  [Ref.  1 and  2]. 
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Least  Squares  Fit  Method  Example 


Given  the  following  set  of  data: 

X 0 1 2 4 7 

f(X)  = Y 0 1 3 12  20 

fit  a curve  of  the  form 

f (X)  = Y = A + BX  + CX2 

Step  1:  Substitute  all  pairs  of  data  into  the  form  equation 
yielding  the  fact  that  the  coefficients  (A,  B,  § C)  must 
satisfy  all  the  following: 

0 = A + B (0)  + C (0)  2 

1 = A + B (1)  + C (1)  2 

3 = A + B (2)  + C ( 2 ) 2 

12  = A + B(4)  + C (4) 2 

20  = A + B ( 7 ) + C ( 7 ) 2 

Now  multiply  each  expression  by  its  coefficient  of  A in 
that  expression  and  add  all  equations  yielding 

36  = 5A  + 14B  + 70C 

Now  multiply  each  expression  by  its  coefficient  of  B in 
that  expression  and  add  all  the  expressions  yielding 
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0(A)  + 0(B)  + 0(C) 


1 

= A 

+ 

IB 

+ 

1C 

6 

= 2A 

+ 

4 (B) 

+ 

8 (C) 

48 

= 4A 

+ 

16  (B) 

+ 

64  (C) 

140 

= 7A 

+ 

44  (B) 

+ 

343(C) 

195 

= 14A 

+ 

70(B) 

+ 

416 (C) 

Now  multiply 

each  expression  by 

its  coefficient  of  C 

that  expression 

and 

add  all  the 

expressions  yielding 

0 

= 

0 (A) 

+ 0(B)  + 

0(C) 

1 

= 

1(A) 

+ 1(B)  + 

1 (C) 

12 

= 

4(A) 

+ 8(B)  + 

16  (C) 

192  = 16(A)  + 64(B)  + 256(C) 

980  = 49  (A)  + 343  (B)  + 2401(C) 

1185  = 7 0A  + 416B  + 2674C 


Now  solve  the  following  three  previously  generated  equations 
for  the  coefficients  A,  B,  & C yielding 


36  = 

5A  + 

14B 

+ 70C 

195  = 

14A  + 

70B 

+ 416C 

1185  = 

70A  + 

416B 

+ 2674C 

A = 

-.99, 

B = 

2.6  , 

and 
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Y = - . 99  + 2 . 6X  + . 065 X2 

The  following  chart  and  plot  depict  the  original  data  and 
the  data  obtained  from  the  equation  for  the  fitted  curve: 


Original 

Fitted  Curve 
Polynomial 

X 

Y 

Y 

0 

0 

-.98 

1 

1 

1.67 

2 

3 

4.48 

4 

12 

10.46 

7 

20 

20.41 

Q • £ . D • 
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APPENDIX  D 


this  page  is  best  quality  practicable 
PROM  COPI  FURNISHED  TO  DDC  — — 


COMPUTER  PROGRAM 


1 H im  " r h i proorhm  Miijri.irr.'.  n.  i;  im  i;  off  ioil  miid  velol :ty.  yuu  mill  iii.ed" 

2 print  " i nr  following  inputs  Pot-:  mi;  conr-uim  ion: " 

3 print 

•i  f i • 1 1 1 r ”1'  ruiiiioy  i-i  i nr  e mi. ri  iube” 

5 PM  ill  luiiinv.'  Ill  IP  III  DEC  F” 

o pr  i n r I'.fiv  •;  ur  min  pmjRLi:  dpi  run  on  or  off" 

7 FRUIT  ”•}'  RELi'iT  I Vi:  HUriI‘1*  I TV 
3 P!  HIT  PUllllil  P III  nil  I tic “ 

,o  pi- hit  ii  IHU  Tt  iri.'i.  r 1 1 in  m lip.;  nun  - 

10  prim  -7  ' MIND  Vl'l  or  1 10  III  KTS’ 

11  IF'  (I  IT  "O'  RUMI-IOY  ■I.OPE  III 
1.':  PMIIT  "0>  Fl  C tO  HI  \ HOC” 

13  PP1II)  “10>r.E.  Fl.FlF'3  SE  T T I MG  III  PEG" 
t-1  PRINT 

15  PRINT  "IE  YOIJ  IIFi'-'E  Fill  THESE  PROCEED . UNFIT  COMPLETE  CIIEOI  YOUR  I ill:  ITER " 

10  PRINT  “fiGlT  1 NOT  Fl.'E  IIOTOPS." 

17  PRINT 

13  DI  3P  "PRESS  OLT  = " j 
l?  INPUT  0 

30  IF  0,3000  THEN  37 

31  IE  00  THEN  37 

33  TiISP  " T UNITlT V TEMP  - ? DEC  F"! 

33  INPUT  BO 
. 4 P — ' 0 0 — 3 3 1 * 0 / 0 
3:.  IP  DO- 130  THEN  13 
30  IF  BO  0 THEN  43 


3 0 

J 1 


40 
4 1 

43 

4 4 
•10 
•1 0 
•17 

40 

Sli 

01 


04 


X::4 1 3 . 030 3000337-0 . 000 1711  300  '-0-3 . 00 04 00 E- 07  •- Fl  1 3-1-3 . 030 1 30  I 
1:0-64-3.4  1057 14E -lO-OIT 

DO  s-n.  0400  34034-7’.  30  3 i 0 1 20  32-00 '0  + 3.70451  1 07E  -00-0  M3 

D 1 •- E- 0 -0.7 03 730  - : 3 • 0 r 3 < •-' . 0 0 >■. 034 0 E - 1 7 + i*l ! 4 

B0=-0.  U01310003  -0 :;-3. 3000 3 1 4 ';F-  07*h+4 . 0733103733E-  lO  oT  +3 

R2-B0-S . 04370077  3C-  1 4 * f ? 3 .-0 . 400  303  SE-  1 8*0  1 4 

B7--1 . 0007E-00+  1 . jOOOi  003 330E-O0*rt-3. 4004  3703E-13-0  f 3 + 3. 04 

B 3=137 - i . 40 1 7137030  IE-1 0*R  + 4 

C -BO+B 1 -L:  i EM *B  T 2 +b  3 *B  I 3 

CO TO  47 

PRINT 

PRINT  "RUNIIOY  OLT  MUST  BE  BTlIII  0 NlID  3000  FT-  TRY  OCR  IN" 
PRINT 

PRINT  ” " 

GOTO  1'3 
PRINT 

PRINT  "TEMP  MUST  EE  DTMN  0 OND  120  DEC  F-  TRY  FiCFlIN" 

PRINT 

PRINT  " “ 

GOTO  23 
PRINT 

PRINT  "F'UNNOY  PRESS  OLT  ="-Fl."FT" 

PRINT 

PRINT  "PUIUIRY  TEMP  =“.BO»"DEG  F" 

PMIIT 

PRINT  — - " 

DI OP  "IF  HD  Oil  FITTER  1-  IE  NOT  ENTER  0"? 

INPUT  D 

IF  D - 1 THEN  00 
IF  " THEN  04 


OE-1  1*013 


337 OE- 1 


Hi  3 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

FROM  COPY  PURNISHED  TODDC  - 


57  PRINT  "PLEASE  ENTER  1 OP  0*  TRY  AGltlN" 

53  PRINT 

5 ■ < r.nro  *.o 

60  E-- 1 . \+  3 780434  + 0 . 50590325*0 

61  pr-im  ■■  double  jin  run  oir 

62  PR  I M l 
£3  GOTO  65 

64  r=0. 54 1777574 0.  55075515*0 

65  PRINT 

66  PRINT  " “ 

67  D I OP  “GROSS  HT-" 5 

68  INPUT  G 

69  IF  G >4 2000  TURN  SI 

70  IF  G< 20000  TURN  SI  • 

71  G 2 = - 4 . 8 S 9 6 E + 0 5 *-8. 4974E+01  *G~5.  7S56E-03*Gt2+l . 9373E-07+G  ? 3-3.  1 744E-12*G»4 

7 2 G 0 -•  G 2+2.0446E-17»  G + 5 

73  G3--5 . 862 1 E+04 - 1 . 01 46E+01 ■ G + 6. SS07E-04^G  1 2-2. 292E-0S*G ' 3+3. 7337E-13+G t 4 

74  G 1 -C 3-2 . 3964  E - 1 3*G  i 5 

75  H-G0+G1+E 

76  PRINT 

77  PRINT  “GROSS  NT  =“iG»"LBS“ 

73  PRINT 

79  PRINT  ” " 

SO  GOTO  86 
81  PRINT 

32  PRINT  “GROSS  NT  MUST  BE  BT1IN  42000  AND  20000. TRY  AGAIN" 

33  PRINT- 

84  PRINT  “ " 

35  GOTO  67 

36  BISP  "FEL  HUMIDITY  = 7."  5 

37  INPUT  I 

33  IF  ICO  THEN  99 
39  IF  I >100  THEN  99 

30  J=0 

31  IF  I <40  THEN  93 

32  J=< 1-40)71000 
?3  K=H* J+H 

?4  PRINT 

?5  PRINT  "PEL  HUMIDITY  =">  l," 

?6  PRINT 

i?  PRINT  " “ 

43  GOTO  104 
>9  PRINT 

100  PRINT  "PEL  HUMIDITY  MUST  BE  BTNII  0 FIND  100.  TRY  AGAIN" 

: 0 1 PRINT 

02  PRINT  " “ 

03  GOTO  36 

04  DISP  “ENTER  RUNWAY  HEADING"! 

05  INPUT  Z 

06  IF  Z/ 360  THEN  135 

07  IF  Z<0  THEN  135 

03  DISP  •ENTER  MIND  DIRECTION  IN  DEG  MAC”! 

09  INPUT  Y 

10  IF  Y/360  THEN  135 

11  IF  Y< 0 THEN  135  • 

12  DISP  “ENTER  WIND  VELOCITY"! 

13  INPUT 

14  W-ABSCY-Z) 
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lie  l - : : • cos <2* pi* 1 1 'ooe > — — ■ - 

llC  II  l.>40  Till.  1 1 1 30 
117'  IT  L 7 -20  1 1 mil  130 
1 1 S ir  L-0  THEN  13-1 

110  LO-G.  7124E+  on  0.  O500E  i.i  1 + I t-2.  330GK--05+I  + 2-1  . 6254 l:  -00+1  I 3-i 3.  3720E--  14*17  + 4 

120  L 1 0000  -1 .0::';.;r-0.  -I.ij.  175U;-0,'-|  1-2-2.  5327'E-l  l-l  t ;s  j . 1 0717-  1 T.- I t ■) 

121  L2 -1.-17S2C-  01- 16G6C-06-I  + 3.4 274C-00  + lvl 2-2 . 78 1 7E- 1 3 • I . t 3 + 9.  3077C- 10* ION 

122  M *L0  + LI *L+  L2*L'»2 

123  GOTO  125 

124  ll-l 

125  nr- 1 nt 

126  FT- HIT  "HEADWIND  = " . L>  “KTS" 

127  PI-' HIT 

123  PI: HIT  " " 

120  GOTO  140 

130  PRINT 

131  PRINT  " HElTDUIlID  MUST  PE  BTIIH  -20  flllD  +40  KTS>  TRY  AGAIN" 

132  PRINT 

133  PRINT  " " 

134  GOTO  103 

135  PRINT 

136  PRINT  "MUG  HEADINGS  MUST  EE  8 THU  0 AMD  360)  TRY  AGAIN" 

13?  PRINT 

133  PRINT  " t " 

130  GOTO  104 

140  Ii ISP  "R/N  SLOPE  - (+  FOR  UP. -FOR  DIO  "5 

141  INPUT  N 

142  IF  U >2  THEN  153 

143  IF  IK •-£  THEN  153 

144  0 0 = 4 . 5 7 0 4 E + 0 1 + 0 . 3 4 2 0 E - 0 1 * N + 2 . 2 2 6 5 E - 0 5 * N + 2 - 2 . 3 3 3 E - 0 9 * M +3  + 7 . 04  1 E-  1 4*M  +4 

145  01-7. 0472-1  1 . 40 1 4E-02+N+0 . 0708E-06  • IT  l 2-7 . 1 235E-  1 0~N  +3+3 . 06S4E-  1 4 «M  l 4 

146  02=5. 3616-8. 5136E-03+M+3.  50 1 4E-0G+N +2-4  . 5032E- 1 U + Nt3+ 1 . 0880E- 1 4*IT+  4 

147  0=00+01 *H+02*N+2 
143  PRINT 

140  PRINT  "RUNWAY  SLOPE 
150  PRINT 

} 151  PRINT  " “ 

: 152  GOTO  158 
153  PRINT 

! 154  PRINT  "R/N  SLOPE  MUST  EE  BTWN  -2  AND  2.  TRY  AGAIN" 

I 155  PRINT 

! 156  PRINT  ” " 

j 157  GOTO  140 

i 153  DISP  "ENTER  A/C  CG<20  TO  35  ‘/.MAO “5 
i 150  INPUT  P 
I 160  IF  P>35  THEN  171 
I 161  IF  P < 20  THEN  171 

1 62  00  = 2. 604 2E  + 03-2.  1 604+0+ 1 . 00 1 5E-03  + 0 + 2- 1 . Ill  OF -07  + 0 +3  + 3. 662E- 1 2+0  + 4 

} 163  01 =-l . 7573E + 02  + 2. 260 1 E-0 1+0-7. 5225E- 05+0 +2  + 7. 70 1 3E-00+ 0+3-2. 5437E- 1 3+0+4 
: 164  02=2. 3540-4 . 0 1 02E-03+0+ 1 . 2332 E- 06+0 +2- 1 . 3234E-1 0+0 +3+4 . 300SE-1 5*0+4 
■ 165  0=00+ 0 1 +P+02+P i 2 
! 166  PRINT 

i 167  PRINT  "CG  =".P.  “’/.MAC" 

163  PRINT 

160  PRINT  ” " 

170  GOTO  176 

171  PRINT 

172  PRINT  “CG  MUST  FALL  BTWN  20  AND  SSMMAC.  TRY  AGAIN" 

173  PRINT  ‘ 

174  PRINT  "— " 

175  GOTO  158 
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rais  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  furnished  to  dlc 


\?e  rn:.p  "enter  te  eu'ip  settiiukso  to  .»o>"» 

. 177  INPUT  R 

: ire.  IF  RMO  THEM  201 
• 1 79  IF  k .20  THEM  201 

| 1 SO  SO 1 . 0079E+02  t 1 . 58ft  1 0254 E -04  v0  12+2.411 1 E-OS  + O t *-:-  S . G737C- 

; 101  01  = 1.  6 1 96E+0 1-2.  -too  51; - 02+0+9.  O'-!  ’ . IF- 00 . 0 ! 2- 1 . 2594E-09  *C'i  2 + 4. 752 

! 1S2  S2-- 1 . 4 758E--01  +2. 359E-04+0-1 . 037C-  07' 01  2+1 . 6016E-  1 1 - 0 1 2-0.310*: 

1 0 3 S = 0 ».h01-R+  S 2 *: R 1 2 
1S4  T 1 = ' im  ■ 100  ' ;•  + 100 

: ISO  T2=S-T1 
■ 186  IF  12,0  THEM  139 

IS?  T=S 

188  GOTO  100 

189  T=<<T1  ! 00 > + 1 ) + 1 00 

190  PRINT 

191  PRINT  "TE  FLOP  SETTING  =" > R> “DEG" 


192  PRINT 

i 192:  PRINT  " " 

194  PRINT  " “ 

195  PRINT  * 

190  PRINT  "ft DJUSTEB  GROUND  ROLL  = " > T i "FT (ROUNDED  UP  TO  NF OREST  10t 
197  PRINT 

193  PRINT  " “ 

j 199  PRINT  " " 

200  GOTO  206 

201  PRINT 

202  PRINT  "TE  FLOP  SETTING  MUST  PE  PTMK  20  AND  40  DEGj  TRY  AGAIN " 

203  PRINT 

204  PRINT  " " 

205  GOTO  176 

206  Ul-5. 4023E  + 0 1 +3 . 4787E-03+G- 1 . 94  75E-03+G  1 2 

207  U=U1  +<26-P;-/2 

208  VG=- 1.9171  E+03+6 . 1 604E+0 1 *11-7. 034SE-01  *IJ  1 2 + 3. 5661E-03+UT  3-6.  6* 

209  VI =7. 6324E+01 -2. 45 1 7-0+2. 8779E-02+U12- 1 . 4753E-C4+U t 3+2. 78?2E-< 

2 1 0 V2  = -7. 2239E-0 1+2. 34 1 5E-02  + U-2. 793E-04+IJ » 2+ 1 . 4596E-06+U1 3-2 . SOI 

2 1 1 V3=V0+V 1 +R+V2+R 1 2 
I 212  V4=  1 NT <V'3) 

213  IF  V3-V4=0  THEN  216 
I 214  V=V4+ 1 
1 215  GOTO  217 

216  V = V 3 

217  PRINT 

218  PRINT  "TAKEOFF  VELOCITY  =“ * V> "KIA3  < ROUNDED  UP  TO  THE  NEAREST 

219  PRINT 

220  PRINT  " *' 

221  GOTO  IS 

222  END 


■13*  0 1 4 
?2E-1  4 -0  r 4 
IE  - 16*0  r 4 


FT) 


57SE-06+U 14 
17*1.114 
’E-09  + UM 


1 KT>- 


LIST  OF  REFERENCES 


A 


1.  Gerald,  C.F.,  Applied  Numerical  Analysis,  p.  284-290, 

Addison-Wesley , 1973. 

2.  Wylie,  C.R.,  Advanced  Engineering  Mathematics,  4th  ed. , 
p.  153-170,  McGraw-Hill,  1975. 

3.  Toumela,  C.H.,  Captain,  USN(RET).  Interview  conducted 
at  The  Naval  Postgraduate  School,  August  6,  1977. 

4.  Winters,  K.A.,  Lieutenant,  USN.  Interview  conducted 
at  the  Naval  Postgraduate  School,  September  30,  1977. 

5 . A7E  Naval  Air  Training  Operational  Procedures  \ 

Standardization  Manual,  NAVAIR  01-45AEE  -1,  15  March, 

1975. 

6 . Hewlett-Packard  9830A  Calculator  Instruction  Manual,  ; 

Hewlett-Packard,  1972. 


51 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 
2 


1.  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Department  Chairman,  Code  67  1 

Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  L.V.  Schmidt,  Code  30  5 

Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93940 

5.  LCDR  W.M.  Siegel  4 

578  Tamarack  Ln 

Lemoore,  California  93245 

6.  Systems  Exploration  Inc.  1 

419  Webster  St. 

Monterey,  Ca.  93940 


52 


